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SUMMARY 

Consideration of the various factors which determine the performance of 
supersonic aircraft and missiles utilizing air-breathing power plants has 
clearly indicated a critical necessity for optimizing the propulsion system- 
airframe relationship. Studies of this relationship demonstrate the dominant 
effects of the required intake and exhaust areas in establishing the airframe 
configuration for high Mach number cruising flight. In the flight Mach number 
range from 2. 5 to 5 0, it is mandatory to consider the power plant as an inte- 
gral part of the aircraft from the standpoint of vehicle lift and drag and balance 
and, in addition, to give careful attention to preserving the maximum possible 
intake efficiency. 

This proposal deals with an integrated ramjet - lifting surface combination 
which is intended to optimize the power plant - airframe relationship mentioned 
above. The essential features of the system are: 

(a) A variable configuration diffuser 

(b) A constant diffuser throat area and Mach number 

(c) A constant burner inlet Mach number 

(d) A variable exit nozzle 

The preliminary analysis contained demonstrates theoretically that the proposed 
system is feasible. The particular system proposed offers the dual advantages 
of high design point performance while permitting efficient off-design operation 
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over a wide Mach number range. The preliminary analysis establishes a 
number of significant advantages characteristic of this system in addition to 
the foregoing benefits. The more important of these include: 

(a) Increased thrust to drag ratio per poimd of payload (particularly 

at high flight Mach numbers). 

(b) Decreased thrust specific fuel consumption for a specific mission 

and payload. 

(c) Structural gains in the form of increased maneuver load capabilities 

and improved flutter characteristics with no associated weight 
penalty. 
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AN INTEGRATED PROPULSION-LIFTING SURFACE SYSTEM 



FOR HIGH MACH NUMBER OPERATION 

PROBLEM BACKGROUND AND GENERAL APPROACH 

The ever increasing demands for speed have forced a three-way compro- 
mise between a thin low drag wing, a wing capable of withstanding high maneuver 
loading, and a light wing easy to design and manufacture. This has necessitated 
increased power and increased demands on target detection and guidance equip- 
ment since the course must be planned at greater distances to avoid high 
maneuver loads. This in turn spells greater weight and a vicious cycle is begun 
with increased cost, materials, and manufacturing skill required at every turn. 
Integration of wing and engine design seemed to offer the best design approach 
for solution of this problem. 

This integration of wing and engine design represents a departure from 
the usual "glider plus power plant" method of design; but the reasoning of L. 

F. Nicholson^* and others seems to indicate that closer integration is a natural 
design evolution as higher and higher speeds are considered. This reasoning 
is paraphrased below; 

Considering the speed range from Mach 0.8 to 4. 0, the most efficient 
aspect ratio decreases as speed increases. The area required for the air intake, 
on the other hand, increases quite rapidly. Considering the frontal area of the 
Mach 0o8 airplane with its large wing span, it was quite reasonable to place the 

* 

The superscripts refer to the authors as numbered in "References". 
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relatively small air intakes and engines rather arbitrarily on a "glider" design. 
In the Mach 4. 0 airplane, however, the situation is reversed and it is quite 
possible that the best design approach might be to "hide" as much as possible 
of the body and/or relatively short wing span behind the air intake. 

In this design proposal the attempt has been made to "hide" the wing 
behind the intake. This was done so the wing could be thickened and thus 
strengthened without adding high drag losses. In the Mach 3. 0 to 5.0 flight 
speed range being considered the choice of a ram jet engine was indicated from 
the operating characteristics of the various types of engines. "Splitting" or 
ducting a thin airfoil and placing a ram jet engine in the wing offered many 
advantages. If the hypothesis is accepted that the most efficient method of 
providing a strong wing is to thicken the cross-section; then the question of 
what should be placed inside the thickened wing, fuel or engine, arises. Wing 
fuel storage is not efficient from the frontal area viewpoint because of the 
rectangular area cross-section and the need for baffles. The wing structural 
members tend to require compartmented tanks which increase plumbing, 
gadgetry, and the possibility of unbalanced wing loading due to trapped fuel. 

On the other hand wrapping fuel tanks around the engine in the fuselage meets 
with combat objections. 

This thickening of the wing would, of course, greatly increase drag. 
However, this drag increase can be avoided in large measure, and further 
weight and drag reduction realized if the engine is removed from the fuselage 
and placed in the wing. Since the engine Intake duct would occupy a portion 
of the frontal area, depending on the flight speed being considered, the small 
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slopes and corresponding low drag of the original thin airfoil are largely preserved. 
Air intake flow also should be considered. The intake ducts should be no longer 
than the length required to efficiently diffuse the air, and the exhaust area should 
be of the same order of magnitude as the intake area. * An engine in the wing 
should reduce this intake and exhaust ducting to the minimum. 

In the NACA publication prepared for the 1957 Triennial Inspection of the 
Lewis Flight Propulsion Laboratory there is an excellent presentation of the drag 
versus flight Mach number of a hypothetical airplane compared to the thrust 
available characteristics of the jet engine. The results of the general analysis 
demonstrate that fixed engine geometry thrust is inadequate, a fixed exhaust 
nozzle is inefficient, and a variable inlet is required to avoid adverse drag 
effects from a spillage pattern; for the operation in the Mach range 0 to 4.0. 

The charts presented in this publication are reproduced in general detail below: 
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If the integrated wing and engine is to operate efficiently throughout the Mach 
range, then, a variable inlet and exhaust is required. 

A cylindrical geometry offers increased total pressure recovery for a 
given free stream Mach number, but the mechanics of providing a variable 



ref, 1 brings this out particularly well in the author's opinion, and also contains 
a good discussion by other members of the Royal Aeronautical Society. 
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geometry while maintaining this increased efficiency lead to overlapping plates 
and "eyelids”; and the whole system is difficult to design, manufacture, and 
maintain. Variability in two dimensional flow, on the other hand, can be provided 
quite simply with light materials easy to manufacture and maintain. In addition, 
the mechanics of a variable cylindrical configuration might well compromise the 
higher theoretical pressure recovery by small interferences to flow, over over- 
lapping plates for example, and thus make the actual pressure recoveries of 
the two systems about equal. If the engine is allowed to extend throughout the 
span of the wing, an ideal means is provided for accomplishing the required 
area variations. 

This study considers the integration of the components indicated above. 

The feasibility of such an integration and the inherent low drag, fuel consumption 
and weight will be demonstrated by considering the performance of the individual 
components; and by calculating the performance of a specific configuration, the 
schematic of which is shown on the next page. The ability of the system to 
withstand high maneuver load while affording a drastic decrease in weight will 
be demonstrated by comparison to a more conventional design. 
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PROPOSED SYSTEM 

Propulsion System 

The basic features of the proposed system which offer significant per- 
formance advantages are: 

(a) a constant inlet throat area, with choked flow; 

(b) a constant burner inlet Mach number; 

(c) a variable area exhaust nozzle. 

This seemed logical from the design objectives since the subsonic diffuser and 
burner efficiency could thus be maximized and the burner design simplified , 

The system proposed employs an Oswatitsch type diffuser to slow the flow 
partially, an "isentropic" diffuser to reduce the flow to Mach 1. 0, the previously 
mentioned constant throat and burner, and a choked flow variable exit. These 
will be considered in order. 

The Supersonic Diffuser 

In the first part of the supersonic diffuser a shock wave pattern is formed 
by the two wedges that make up the shock generator; and external compression 
takes place, followed by "isentropic" internal compression to Mach 1,0. 

Since only two large shocks are used in the external compressor, there is a 
loss in efficiency over the ideal completely "isentropic" diffuser. However, 
by accepting this initial loss, three advantages are gained. They are: 

(1) A smaller initial capture area is required to provide full flow in the 
throat. This occurs because of the following: Assume that a pure "isentropic” 
diffuser were used throughout the supersonic flow. Since losses would be 
reduced, the capture area would have to be increased to provide full flow for 
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a given size throat. If the burner cross-section is to remain the same, this 
would require large flow expansion angles on the outside of the diffuser lips at 
high speed and thus greatly increase the drag. If the burner cross-section is 
increased to reduce these angles, the drag near Mach 1,0 is greatly increased. 

If on the other hand, the duct required is shortened to a minimum by using only 
one or two oblique shocks and a finite strength normal shock, the external 
compression limits would be exceeded or else high losses will be suffered 
across the normal shock. In either case compromising losses would be suffered. 

It is Indicated, therefore, that some loss must be accepted; and the shape 
of the supersonic compression curves, as well as comparative losses across 
the normal shock, and space and weight requirements strongly suggest 
acceptance of loss in the high supersonic flow region. 

(2) A shorter duct is required with less boundary layer, drag, weight, 
and complexity. 

(3) An "isentropic" diffuser required a zero angle leading edge It will 
be shown later that the angle required for the leading edge for a design speed of 
Mach 3. 0 is approximately ten degrees. It was felt that this angle could be 
realized in manufacture and would have sufficient strength to retain its shape 

in flight. 

The use of external compression offers a great advantage by reducing the 
throat area necessary for swallowing the shock when starting the engine. This 
normal shock will be present at the forward part of the lip at the entrance to 
the internal supersonic "isentropic" diffuser when the engine is not producing 



thrust, due to the difference in flow conditions between an engine that is producing 
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thrust and a cold engine. When the engine is started this normal shock will be 
swallowed (move aft of the throat) and then be forced up to the throat by the 
action of the variable exit nozzle that will be discussed later. The fact that 
the capture area is variable permits easy starting of the engine, and allows the 
engine to operate at peak efficiency since the ratio of the capture area to the 
constant throat area can be changed from that required for swallowing to that 
required for peak efficiency. A schematic cross-section of the proposed system 
is given on the next page. 

The Subsonic Diffuser 

Ref. 8 gives the relation of total pressure recovery to the rate of change 
of the cross-sectional area per foot of length of a subsonic diffuser. It can be 
seen that this curve peaks quite sharply and has a definite maximum. Since the 
diffuser throat area and burner area are kept constant, the subsonic diffuser 
can be designed to realize this peak efficiency for all values of free stream 
Mach number. This, of course, would not be possible without some difficulty 
if the throat area were allowed to vary. 

The Burner and Exhaust System 

Because of the constant burner Mach number, the burner design can be 
simplified for low cost and ease of maintenance and optimized for peak efficiency; 
and this efficiency will be maintained for all values of free stream Mach number. 

A two dimensional variable exit is proposed (see Fig. 7) to accelerate the 
flow to the local speed of sound at the throat (choked flow) and provide expansion 
of the flow to ambient pressure. This makes the engine easier to control and, 
as will be shown in the calculations, increases the thnist as much as 60%. Because 
of the large thrust to drag ratio, particularly at high speed and high altitude, of 
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the proposed wing engine combination; cruise will be at reduced throttle. The 
throat will have to close further under this condition to maintain choked flow. 

At some point, then, the variable throat will be positioned so that the exit area to 
throat area ratio is optimum and ideal expansion will result with greatly increased 
thrusto The improvement in performance that is present over the entire speed 
range will be greatest at this point. 

Lifting Surface 

It will be shown later that for best total pressure recovery when operating 
below the design speed, the shock generator must be rotated up. This rotation 
places the outer surface of the shock generator at some angle to the relative 
wind. It was found that the lift produced from this free stream deflection could 
provide the lift required for level flight at moderate altitudes while the total 
pressure recovery of the inlet air was kept at or near the maximum. This is 
not the best way to produce lift, dragwise, since it involves larger angles over 
shorter chord lengths. It is, however, a very efficient method when air intake 
and compression is also considered, because of the double use of the shock 
generator. The relative efficiency of ”by passing” some of the air not required 
for engine operation at below design speed, converting it to lift, and accepting the 
relatively high drag versus allowing the excess air to enter the intake and then 
trying to dump it with as little loss as possible will be discussed later. Structur- 
al benefits that also accrue from this action will be discussed in the next 
section. 

Structure 

For structural purposes the usual supersonic wing is only considered 
in the area between 25% and 75% of the chord since the pointed leading and 
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trailing edges outside of this area provide so little resistance to bending that 
their vertical load supporting capabilities are negligible. With the system 
proposed, however, the effective moment arm is available over an increased 
portion of the chord. This provides an increased load carrying capability of 
the wing that is further increased because the full maximum thickness of the 
normal wing can not be used as a moment arm, even between 25% and 75%, 
because of upper and lower wing surface cuivature. 

Even a greater margin of strength is provided, however, by the large 
increase in wing thickness and thus moment arm. The recent design problems 
and special manufacturing techniques are thus avoided, and wing design is 
reduced to standard techniques where a great deal of background data is available. 
Even considering the additional weight of the engine inside the wing, this wing is 
easier to design and manufacture, lighter, and capable of supporting more 
vertical load than the normal thin supersonic wing. This is particularly 
advantageous for high altitude intercept since it reduces the requirements on 
intercept angle, flight path, detection equipment, and computer equipment; 
and takes the latter two items out of the long range radar and crystal ball 
category. Due to the lessened requirement or intercept angle, early wai'ning 
detection requirements are reduced since an intercept is possible by flying in 
the direction of the target instead of going around behind and overtaking. This 
shortens the flight path and reduces time and weight requirement. 

Aeroelastic effects that cause the wing to twist when an aeleron or other 
control surface is deflected are greatly reduced because of the enlarged "torque 
box" (increased wing thickness and usable chord length) that is available. The 
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elastic deflection when high maneuver load is applied is greatly reduced as 
previously mentioned. Therefore aeroelastic effects weren't considered. 

The effects of inertia coupling were not computed, but it can be seen that 
the thickening of the wing and the reduction in size and weight of the fuselage 
will give a much more favorable mass distribution. The proposed system has 
been p>artitioned by solid ribs in order to support the skin and stringers, elimi- 
nate any transverse wave formation, and, by extending these ribs forward, to 
provide a secondary means of supporting and positioning the shock generator 
and compressor lip. This is shown on the oblique view previously given in 
the "Problem Background" section. 

The structural load on the compressor lip of the proposed system is not 
as great as might be imagined. The outside of the lip is inclined to the relative 
wind; therefore, there is a pressure force behind the shock wave it creates that 
tends to push the lip in. The intake flow that has been compressed by the shock 
generator tends to push the lip out. At the leading edge of the compressor lip 
these two forces will nearly balance each other and result in a small moment 
even though the moment arm is relatively long. As the compressor lip turns 
the flow, the pressure pushing the lip out will begin to dominate the force 
balance; but at the same time the moment arm will be decreased and thus a 
relatively low moment will be maintained. 

If this pressure balance were ignored, one might well attempt to decrease 
the drag of the system by somehow making the outer surface of the compressor 
lip parallel to the free stream. However, from the above it can be seen that 
this would force a drastic increase in moment on a member that must be kept 
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thin for best efficiency. Therefore, the net result of this attempt might well be 
a loss in net thrust (thrust minus drag) through loss in compression efficiency 
as well as an increase in weight. 

This balance of pressure forces also works to advantage in the case of the 
shock generator when it is positioned for operation below the design speed. 

Here again tiie outer surface is inclined to the relative wind and the inner and 
outer pressure forces tend to balance each other. 
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SYMBOLS 

A area (feet^) 

b wing span (feet) 

C Dimensionless coefficient, also wing structural chord 
c wing chord (feet) 

Cj’ thrust coefficient based on dynamic pressure 

Cf exhaust thrust coefficient 

D total drag (pounds) 

d the couple moment arm to resist wing bending 
dj^ hydraulic diameter (feet) 

F thrust (pounds) 

g acceleration due to gravity (32. 2 feet/ sec ) 

K thousands (lOK is ten thousand) 

L total lift (pounds) 

M Mach number 

p pressure (pounds/feet^) 

total pressure recovery ratio of the supersonic diffuser 

R body radius (feet) 

2 

S wing area (feet ) 

S. L. Sea Level 

T temperature (°R) degrees Rankine 

V velocity (ft/sec) 

VV airplane weight (pounds) 

\v engine mass flow (pounds/sec) 
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Greek leters 

angle of attack (degrees) 

angle of the shock wave relative to the flow in front of the shock (degrees) 
specific heat ratio 
efficiency factor 

wedge angle of flow deflection (relative to the flow in front of the wedge) 
(degrees) 

angle of deflection of Prandtl-Meyer "flow around the corner" (degrees) 
exhaust exit area versus throat area ratio 
Subscripts 

0 free stream conditions 

1 conditions after the first shock in the diffuser when considering the diffuser 
only, and at the diffuser throat when considering the complete engine. 

2 conditions after the second shock in the diffuser, and at the entrance to 
the burner when considering the complete engine 

3 conditions after the third shock wave, and at the exit to the burner 

D drag 

DP design point of the diffuser 

e exit conditions at the exhaust nozzle 
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y/ 



V 



L 

N 

T 

W 



lift 

Normal shock 

stagnation or total pressure (^Tq stagnation free stream pressure) 

(pounds/feet ) also the exhaust nozzle throat conditions, 
wedge angle relative to free stream intake flow when the shock generator 
is positioned for some Mach number other than (degrees) 
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ANALYSIS 

The analysis of this proposal has been divided into two parts; the outside 
of the wing and engine system, to include aerodynamic and structural con- 
siderations; and the inside, or ducted area, to include the supersonic diffuser 
scheme and engine considerations The thickness of the overall wing was fixed 
by the required capture area of the air intake which is in turn dictated by the 
throat heighto The chord of the wing was fixed by the following factors: 

(a) the required length of the shock generator or supersonic diffuser 

(b) the length required by the subsonic diffuser so as not to exceed 
efficient diffusing rates 

(c) the length required for the burner 

(d) the length required for the variable exit 

By expressing the burner length and exit length in the customary burner 
diameters (in this case burner height? it can be seen that the chord length is 
also dictated by the throat height, if the burner Mach number and thus area is 
kept constant. 

In this section only the drag and thrust per span foot of wing will be con- 
sidered. It will be further assumed that the wing has a constant cross-section 
A diffuser throat height of two inches was arbitrarily chosen. Since all wing 
cross-section dimensions are based on this throat height, all dimensions were 
specified. This made the actual chord: 

Shock generator 1 625 ft. .271 Total chord 

Aft wing 4,375 ft. . 725 Total chord 
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The height of the throat is limited by the length of the generator that can be 
controlled and held in place, but 1 625 feet seemed very reasonable. Using 
the optimum diffuser configuration that will be determined for M^p 3.0 the 
performance of one spanwise foot of wing was determined by the calculations 
as shown in Appendix I and II. The method of calculation used and the assumptions 
made will be discussed in order. 

Determination of Optimum Supersonic Diffuser Configuration 

In the original theoretical analysis of this type of diffuser by K Oswatitsch^ 
it was found that the maximum total pressure recovery through n-1 oblique shocks 
and one normal shock is accomplished by making the total pressure ratios across 
all of the oblique shocks and the normal shock equal This criteria will be used 
to find the strength of the two oblique shocks for optimum performance of the 
first portion of the supersonic diffuser. The second portion of the supersonic 
diffuser employs "isentropic" compression. "Isentropic” has been placed in 
quotation marks thus far to denote that there are some losses in the compressor 
and thus the flow is not actually reversible. These losses will be discussed 
later and shown to be small, thus the flow does approach the isentropic condition. 
This flow is actually Prandtl-Meyer flow around a comer or zero radius turning. 
Parameters for this flow can be found in any good aerodynamics book.^ 

The flow through the diffuser throat is parallel to the free stream flow. 
Therefore, the angle of deflection of the Prandtl-Meyer flow is equal to the total 
angle of deflection of the flow caused by the two oblique shocks (see Fig, 1), 
Prandtl-Meyer flow deflected through some angle has a unique Mach number. 



By considering the flow to be compressive instead of the normal expansive flow, 
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this unique Mach number will be the Mach number before deflection. Parameters 
for this type of flow (diffuser shape, focal point, etc.) are given in TN 3589^ and 
have been confirmed by experimental results Matching this unique Mach number 
of the flow before Prandtl-Meyer deflection to the Mach number of the flow after the 
two oblique shock deflections by iteration, and using the Oswatitsch criteria for 
equal total pressure losses across the two oblique shocks gives a unique optimum 
configuration for diffusing any given free stream Mach number to sonic velocity 
Fig. 2 gives; as a function of the free stream Mach number, M^; this optimum 
configuration in terms of the first wedge angle of deflection to the free 
stream; Mj^, the Mach number of the flow after deflection through degrees; 

0^, the second wedge angle of deflection to the previously deflected flow; M 2 , 
the Mach number of the flow after the second deflection; \J , the total angle of 
oblique shock deflection equal to the deflect! in angle of the Prandtl-Meyer flow; 

the generated shock wave angle relative to free stream M^; and 32, the 
generated shock wave angle relative to the flow after the deflection by 3^. These 

C 

values were taken from charts. By applying the condition that the two shock 
waves intersect at the isentropic diffuser lip, a unique configuration is specified 
by some given throat height. Fig, 1 was constructed, as an example, for a 
optimum diffuser of M^ 3. 0 flow where the throat height was taken as one inch. 

This configuration of the shock generator is, therefore, defined as Mach Design 
Point for a free stream Mach number of 3. 0, M^p 3. 0. Considering only the 
optimum configuration for each free stream Mach number as given in Fig. 2 (Mj^p 
and thus the configuration changing), the total pressure recovery ratio between the 
free stream and the diffuser throat was computed from charts^ to give the maximum 
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theoretical total pressure recovery for the proposed configuration at any given 
Mop- The results are given in Fig, 3 for a two oblique shock system. The 
improved performance of a three shock system will be given later. 

Design speed pressure recovery curves, adjusted for losses as discussed 
in Appendix Ht, have also been included in Fig. 3 for optimum configurations 
using three oblique shocks plus one normal shock, two oblique shocks plus one 
normal shock, and one oblique shock plus one normal shock. Comparison of 
these curves show the proposed design to be superior at design speed, 

The Constant Configuration Shock Generator 

If the optimum configuration for some Mpp is taken, and a constant 
configuration shock generator constructed and rotated about the point at the 
focal point of the isentropic compressor (Point A in Fig, 4), a unique position 
of the shock generator relative to free stream may be determined in the following 
manner: 

At some rotated position of the shock generator the outer surface of the 

shock generator will form some wedge angle relative to the free stream, 

The parameters of this outer or by-passed flow will be discussed later. The 

inner surface of the shock generator will also form some wedge angle relative 

to the free stream intake flow, ^ . Since the shock generator configuration 

w 

is constant: 

Q *9 = 0^ z constant 

o w 1 

n 

^ z constant 
2 




For some given is specified and, since B 2 is constant, M 2 and V are 
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also specified. By iteration to match M 2 and V to the Prandtl-Meyer flow 
characteristic Mach number and angle of deflection, a unique position of the 
shock generator can be determined for any given free stream Mach number. 

This of coarse requires that the isentropic compressor be retractable, since 
a different deflection angle is required for each M^. The isentropic compressor 
lip, then, is a preformed member that is retracted into a channel as controlled 
by the same free stream Mach number sensing device that controls the position 
of the shock generator. The isentropic compressor profiles given in ref, 4 
demonstrate that this constant shape of the lip is possible over the range of 
M 2 being considered (M 2 *= 2,22). The advantages of this movable lip have 
been mentioned previously and will be discussed in more detail later, as will 
the structural considerations. This scheme of operation is shown schematically 
in Fig, 4 where the position of the leading edge of the lip is noted by the niunber 
corresponding to the free stream Mach number, M^. 

Operation of this configuration below design speed will result in slightly 

less than the given ^T^/p pressure recoveries, To demonstrate the 

max 

magnitude of this loss, the same design speed originally taken of Mj)p 3. 0 was 
again selected. The theoretical pressure recovery curve of this configuration 
is shown in Fig, 5 and labeled "Constant 0 ^ and for Mj^p 3. 0. It can be 
seen that this performance is slightly less than the maximum possible for Mach 
numbers less than 2. 0 of other configurations when they are operating at design 
speed, but is suj>erior from there to beyond design speed The AIA Ram Recovery 
Curve has also been included for comparison. Fig. 5 shows the loss in total 

p 

pressure recovery as compared to T^/p 

'^max 



for operation at speeds below 
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The partial external compression and partial internal compression 
provided avoids the flow deflection problem at the exterior of the isentropic 
compressor lip and allows the lip to have a finite angle at the leading edge 
without causing a detached bow wave and the accompanying spillage and drag 
This can best be seen by noting that slightly subcritical flow is easier to control 
and increases pressure recovery characteristics slightly Therefore, best 
operating conditions would move the shock waves caused by the shock generator 
slightly forward of the critical flow condition shown in the figures. With this 
shock pattern moved forward it can be seen that the isentropic compressor 
lip lies entirely within the deflected flow caused by the shock generator. The 
leading edge angle of this lip can, therefore, be made some small angle without 
producing a detached bow wave, 

Boimdary Layer Considerations 

As shown in Fig. 6, it is assumed that the bovmdary layer will begin to 
build on the inner surface of the shock generator This layer will tend to pull 
away from the surface as it nears the flow pattern caused by the isentropic 
compressor lip. The ducting at A is such that it scrapes off a portion of that 
layer before the pull away tendency becomes pronounced. The boundary layer 
scraped off is utilized for cooling as described later. The same scheme of 
scraping off the boimdary la}'^er on the outer surface of the shock generator is 
used at Point B (see Fig. 4). Wind tunnel experimentation has shown that the 
boundary layer on the isentropic compressor lip can be minimized by holes 
drilled at an angle in the compressor lip. These holes are positioned so that 
the number exposed to the flow will be increased in proportion to the tendency 
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for the boundary layer to grow as the compressor lip is extended. 

Subsonic Flow 

The flow thus enters the throat [see Figo 4) with a very small boundary 
layer at a speed very close to sonic velocity. After the engine is started, a 
normal shock wave is formed at the throat and kept there by burner and exit 
conditions that are controlled by the variable exit that will be discussed later» 
Therefore, choked flow with constant diffuser throat area is maintained at all 
engine operating conditions. Since the flow change across the normal shock is 
from very slightly supersonic to very slightly subsonic, a minimum total pressure 
loss is experienced across this normal wave. The subsonic diffuser losses can 
be minimized by the use of the optimum expansion angle in the diffuser. Space 
has been provided to efficiently diffuse the flow to Mach 0. 2 at the burner, and 
a conservative efficiency of . 90 was assumed for this operation in the later 
calculations. 

The Burner and Exhaust System 

The burner has only been shown schematically. As mentioned previously, 
a simplified design is possible because the inlet Mach number is maintained 
constant. In thrust computations it was found to be convenient to use a 
dimensionless coefficient, Cp-, that was computed by: 

C = T 

^ ^2 Pq M§ C 

where C is the average wing chord and T is the thrust of one lineal foot of wing 
These computations were also made in Appendix I and II, and the computed values 
of Cp are given on Figs, 10 and 14. 



23 . 



Since exit axea versus throat area ratio, C, was used to some extent in 
the calculations, the usual ram jet thrust equation was rewritten as: 

I' = ^e - 

where Fg = Cj Pj Aj and Fq = Vq w/g 

3 

7 

values of were then obtained from charts. The area of the exit is limited 

by the dimensions of the wing. When the optimum t was exceeded (about 2. 2) 

the area of the throat was decreased to make this ratio optimum. This in turn 

decreased the fuel air ratio for operation well below stoichiometric. 

It was found from a determination of Cp that at low supersonic 

speeds C rises with altitude, but this rise is limited by the total temperature 
r 

rise at stoichiometric fuel ratio (thermal choking) above some optimum altitude. 

It was assumed that 20, 000 feet altitude provided the maximum Cp/ Cjj ratio. 

At this altitude the shock generator produces a very large portion of the lift 
required for level flight as previously mentionedo Only the altitude of 20, 000 
feet was considered in the low supersonic speed range, since showing thrust to 
be well over the drag for one altitude demonstrates the ability of the vehicle to 
accelerate to cruise speed. 

Cooling 

Cooling is accomplished by the two streams of boundary air bled from 
Points A and B as previously mentioned. The air from Point A is channeled 
through hollow posts that make up the middle portion of the front spar as shown 
in Figs. 4 and 8. This flow is then led aft to cool the compressor lip side of 
the wing, and exhausted out the rear. Air from Point B is used to cool the 
shock generator side of the wing. Increasing the size of the wing would increase 
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the size of these cooling passages accordingly. 

Lifting Surface 

Since the engine air intake has been located in the leading edge of the wing 
and accounts for part of the air-stream ahead of the wing, the drag-wise effect 
of thickening the wing is greatly modified. In order to estimate the magnitude 
of this effect the dimensionless lift coefficient, was computed as the lift 
required for level flight from the formula.- 



where VV/S was taken as a constant 200 pounds per square foot; and the dimension- 
less drag coefficient, C^j, was computed for one foot of wing using a diffuser 
throat height of two inches, which gives a sectional dimension of the shock 
generator per 24 square inches of throat area. Refering to the figure below, 
the curved portions of the upper and lower surfaces (X and Y respectively) were 
approximated by extending the sloped surface to the intersection with a horizontal 
line. Since the diffuser has a practically constant orientation with the relative 
wind, the horizontal portion of the simplified model will have no drag. 
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It was felt that this simplification would yield results certainly within the 
accuracy of the method of linearized thin airfoils used to compute the drag. 
The outer surfaces of a diffuser for a given then, have a chord-wise 

length and slope for both the upper and lower surface that is dependent on 
Linearized thin airfoil theory was then utilized to obtain a Cj: 



CT^ = 



D ■ '-aveiyMg - 1 ^ 



L 2 

+ 

l dx (up 



dx 



dx 



low 



where L is the chord- wise length. It should be noted however, that this is 
based on an average chord length from the tip of the shock generator to the axis 
of rotation of the shock generator. When a of the remaining portion of the 
wing is computed based on its chord, the two chords must be totaled and the two 
Cjj’s combined in proportion to their characteristic chord's fractional part of the 
total chord. Thus if the average chord length of the shock generator for some 
Mdp was 1/4 of the total chord, the total would be; 



1/4 C + 3/4 C. = 

“(shock “(remaining ^otal 

gen. ) wing) 

A based on linearized thin airfoil theory was also computed that must be 
combined in the same manner: 



C 



L 



40 : 




[ 57.3 



(O in degrees) 



The above value must be divided by two when computing the shock generator lift, 

since only one surface is considered. C and C , were computed for the remaining 

L d 

wing by the normal application of linearized thin airfoil theory. The lift of the 
intake air passing through the engine was assumed to be zero by the following 
argument: It will be shown later that the aft portion of the wing remains practi- 
cally constant at zero degrees relative to the free stream. The vertical component 
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of the momentuin flux entering the engine is zero since the free stream is 
horizontal. If the angle of attack of the aft portion of the wing is zero, then the 
vertical component of the momentum flux leaving the engine is zero. Therefore, 
the vertical momentum change and thus the lift must be zero. 

Computations were made in Appendix I for a two shock system with 3. 0, 
and in Appendix n for a three shock system with M^p 4. 0. In some cases angles 
greater than the normal limit of accuracy of the formulas were used. However, 
in each case the proportional part of the total term was considered so that the 
resultant error was reduced to an acceptable value when considered as an error 
for the whole term. The computed values of are given on Figs. 9 and 12. The 
results of these calculations are used as a basis for estimating total drag in 
order to demonstrate the feasibility of the system. 

It can be seen from Fig. 4 that as M is decreased; S , the outer wedge 

o o 

angle of the shock generator, is increased. Therefore, there will be a speed 
at which the shock wave attached to the leading edge of the generator tends to 
become detached and form a bow wave, with the flow behind this wave subsonic. 

At this point, if is decreased slightly, the shock wave will again become 
attached to the leading edge, but the normal shock will move from the throat 
to the lip of the isentropic diffuser, since the turning angle inside the diffuser 
is too great. vVith the normal shock wave positioned at the lip, it can be seen 
that the isentropic diffuser is now acting as a normal shock diffuser with the 
Mach number ahead of the normal shock greater than a slight amount above Mach 
1.0. A total pressure loss will now be incurred across the normal shock at 
the lip. 

During this regime (when the leading edge shock wave is kept attached to 
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the leading edge by decreasing 0 ^ ) and will increase more slowly as 

is decreased. This retarded increase in will force an increase in for the 

remainder of the wing to maintain level flight; and some point will be reached, 

depending on the individual configuration, where Cj. / Cj^ will become a maximum 

for each Mq. The location of this exact optimum point for each M^, using 

theoretical data is rather meaningless, since it involves the difference of relatively 

large quantities to find a minimum. These exact optimum points are not critical 

if it can be demonstrated that at each point the thrust is well above the drag and 

thus the vehicle is capable of acceleration through this region. 

In this investigation the position of the generator was arbitrarily chosen 

as shown in Fig. 9 by the following schedule; 

^ 2. 2 Best total pressure recovery (shock in the throat). 

2. 2 =* M>1. 5 Limit of external compression (shock on the lip). 

o _ 

This compression limit is shown on Fig. 9 in terms of M£ as This value 

was taken from NACA TN 3589. The indicated in the figure was slightly 

increased to insure a stable shock system. The position of the shock generator 

required to produce this shock system is given on Fig. 9 in terms of "9 , M 

o o 

less than 1. 5 was not investigated because theoretical results in this transonic 
region have doubtful accuracy, and are properly a case for wind tunnel investigation. 
Considerations for operation of the compressor at speeds below this region con- 
sidered above will be discussed later. At this point it suffices to say that for the 
schedule chosen, no normal shock will form either on the inner or outer surface 
of the shock generator. Since there will be no shock detachment, no addition 
drag from intake air spillage will be present. 
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In a later section an approximate solution is obtained for the proposed 
system with a body and tail attached in order to compare the proposal with a 
conventional vehicle with a supersonic thin air-foil wing. 
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DISCUSSION 

Operating Characteristics 

In order to demonstrate the feasibility of the proposal, the previously 
calculated values of C^j and Cp. per spanwise foot of wing were used to obtain 
the performance of a vehicle using the proposed system.. These calculations 
were also made in Appendix I and n. The resultant thrust and drag coefficients 
have been plotted versus Mach number for various altitudes in Fig. 10 for the 
two shock system with Mpjp 3. 0 and in Fig. 14 for the three shock system with 
Mdp 4. 0. It can be seen from these plots that the large amount of excess 
thrust provides excellent maneuverability in a wide operating range that will 
exceed the altitudes and speed values computed This large amount of excess 
thrust will dictate reduced power settings in most cases, so that the burner 
temperature calculated will probably not be obtained except while operating 
close to Mach 1. 0. The scheme of maneuver here is to traverse this area as 
rapidly as possible and to avoid this range in the same manner that most aircraft 
avoid certain operating areas. 

The effect of wing moment coefficient was not considered. Wind tunnel 
tests will probably prove that this configuration can be operated through the 
transonic region and down in the subsonic region^ However, the theoretical 
calculations required for sufficient accuracy in these areas are beyond the 
scope of this proposal and the ability of the author. The spanwise length of the 
wing was fixed by consideration of such factors as the wing loading, the aspect 
ratio, and the proportional part of the drag caused by body and tail surfaces 
that each foot of engine would have to account for„ Because of the use of 
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dimensionless coefficients of drag and thrust and the unique relationship of all 
cross-section dimensions, the Cj. and Cj-j values determined hold for any throat 
height as long as the wing loading is held constant. Therefore, if the body weight 
and shape are held constant, the aspect ratio is inversely proportional to the 
throat height and can be set from considerations of best structural efficiency. 

A wing loading of 200 pounds per square foot was chosen. 

The closed body shape given by Haack was used for these calculations 
since it has a minimum drag for a given volume and length: 

Profile: R/L = Rj^^j^/L(x/L - 
The approximate wave drag of the body is given below: ^ 




C per ft. 
D 

of wing 
for a 6 ft. 
chord 



1.5 2.0 2.5 3.0 3.5 4.0 



Assume that the drag of the rest of the body and the skin friction of the complete 
body equals the body wave drag. Assume that R^„^ = 1.38 Ft ; the frontal area 

XII LiX 

p 

will then be 6. 00 ft. , and the equivalent frontal area to include the rest of the 

2 

body and skin friction will be 12. 00 ft. . Referring this to a lineal foot of wing 
with a chord of six feet will multiply the above Cj) value by two. If the lei^th 
of the wing is taken as six feet, each lineal foot of wing will have to account for 
l/6 of the above Cjj or l/3 of the given Cj) value. This scale is plotted on the 
right side of the figure. This gives, on the following page, the values for Mq 



considered: 
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^0 


1.5 


1.75 


2.2 


2.6 


3.0 


3.5 


4. 0 




.0202 


.01975 


. 019 


. 0182 


.0175 


. 0165 


. 0155 



It was felt that this was an overly large body and all assumptions conservativeo 

The actual span of the wing will, therefore, be 8.76 feet and the aspect ratio 1. 46. 

It was assumed that the body produced no lift, the vehicle had a gross weight of 

about 7, 500 pounds (300 lb. fuel assumed consumed accelerating to Mq = 1. 5). 

The Cn and Cv calculations shown on Fig. 10 were made using the 
^ ^ max 

variable exit nozzle previously discussed. It was assumed that the missile 
could be rolled 180° (turned upside down) at cruise speed in order to utilize 
the lift of the isentropic compressor lip at 3. 0, and the lift of the shock 
generator at greater speeds. The probable behavior at speeds between 1. 0 
and 1. 5 has been indicated on the figure in dashed lines, on the basis of typical 
operation in this region. 

The angle of attack was found to remain nearly zero (0° to 6°) for this 
condition. The condition where the vehicle was restricted from rolling was also 
investigated. The angle of attack was found to be greater (1. 2° to 7. 7°) for this case. 
However, in both cases Cp remained well above as can be seen in Fig. 10. 

Loss of lift at the wing-tip was not considered. However, since most of the lift 
is produced at the forward part of the wing, a very small portion of the lift 
producing surface lies within the wing-tip Mach cone. This loss of lift was, 
therefore, considered negligible. 

Performance of a Three Oblique Shock System with M p p 4. 0 

A three oblique shock system was investigated for a Mach Design Point of 
4. 0. The configuration shown on Fig. 11 was taken from the chart on Fig. 12. 
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Fig. 12 also gives the limiting Mach number for external compression as 
discussed previously; and the total pressure recovery of the supersonic dif- 

P 

fuser, » when the diffuser is operating at design speed, MjQp. It 

should be noted that the total pressure recovery is greater than the theoretical 
maximum for a fixed geometry isentropic diffuser with maximum external 
and maximum internal compression. This maximum is shown on Fig. 12 as 
given in ref. 4. Exceeding this curve is possible because of the variable 
configuration of the proposed design. This can be readily seen by recalling 
that the external compressor is limited as previously mentioned; and in the 
theoretical maximum case by the area ratio of the internal compressor entrance 
to throat area for shock swallowing. This limit does not exist for the proposed 
design since this ratio can be controlled by moving the small lip. 

The schedule for the position of the shock generator for this study was 
arbitrarily chosen as shown in Fig. 13. It is summarized below; 

3. 0 Best total pressure recovery 

3.0=* 2. 0 Arbitrary schedule (The curves shown on P'ig. 12 and 

13 could be considerably straightened by selecting the 
optimum shock generator position. ) 

2. 0 > 1 . 75 M 3 slightly greater than Mj_^ (the limit of external 

compression). 

1. 755 -Mq=» 1. 5 M 2 slightly greater than M^. 

It is previously stated that methods to extend the low supersonic range 
operation without a normal shock for ming either on or in front of the shock 
generator; and the relative efficiency of "by passing" some of the air in front 
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of the wing, converting it to lift, and accepting the relatively high drag; would 
be discussed later. These can best be illustrated by the design now being con- 
sidered in the Mach 1.5 to 1. 75 range. Here it is proposed to hinge the last 
ramp and make zero at speeds below 1.75, as shown in Fig. 11. This 
allows Mg to approach M 1. 0 without a normal shock forming at 0 ^. This is 
directly applicable to the two shock case by hinging the last ramp there. This 
positioning device need only be a two position system and can thus be kept small 
and light. One is tempted to utilize this hinged ramp system at greater speeds, 
and thus be able to hold the position of the shock generator fixed, and acliieve 
design simplification. Consider, however, the thrust equation; 

Pe = Po 

we Wo 

F z Vg - (optimum exhaust expansion) 

The V term increases as the square of the velocity. At low supersonic 

O 

speed some a dditional air can be utilized for cooling between the burner wall and 
wing skin; and then exhausted at the trailing edge to recover this lost thrust. If 
the variable ramp Is attempted exclusively, the passages are no longer capable 
of handling the large mass flow; and air must be dumped, with an accompanying 
high loss of thrust, 

Cp. and Cp calculations were made as before in Appendix n and are 
^ -“^max 

given in Fig. 14. The same body shape and dimensions and wing loading were 
used. It was felt that the increased weight enclosed in the body was more 
realistic, but still conservative. The height of the diffuser throat was taken 
as one inch, making the diffuser length 27 inches and the aft wing 27 inches for a 
total chord of 4. 5 feet. The length of the wing was taken as 12 feet, giving a 
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wing span of 14. 76 feet, an aspect ratio of 3. 28, and a gross weight of the 

vehicle about 11, 000 pounds. Again it was assumed that the missile could be 

rolled for cruise. The angle of attack was again found to remain nearly zero 

(0° to 1. 23° at 20, 000 feet) for this condition. The effect of restricting this 

roll was not investigated, since it was shown previously to be negligible. Again 

Cp, remained well above as can be seen in Fig. 14. Wing-tip loss of lift 

was not considered for reasons previously mentioned. For this design the relative 

thickness of the wing increased while the relative thickness of the burner remained 

constant. Increased t ratios were, therefore, possible and a value as high 

as 1. 596 was realized at M 4. 0. 

o 

The basis for determining the span was the requirement of a satisfactory 
margin of excess thrust for vehicle acceleration in the Mq i. 5 region. If partial 
boost can be permitted in this range, a further extension of range and/or de- 
crease in weight is possible. 

It is shown in Appendix I and n that in both 3. 0 and 4. 0 cases, CX 
remains very close to zero throughout the operating range considered until 
cruise speed is reached. If the missile is rolled, CX is again very close to 
zero for both cases in the cruise condition. 

Comparison of the Proposed Design to a Conventional Design 

Since relative wing strength and ability to withstand high maneuver load 
will be an important part of the comparison, the conventional configuration 
selected for comparison should have the engine located on the wing and not in 
the body. This will remove the body and wing root strength from consideration 
and allow attention to be focused on the wing itself. To further aid in the comparison. 
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wing area plan shapes should be the same; that is straight, low aspect ratio 
wings. To further simplify the problem, consider the conventional engine to be 
supported on the wing in such a manner that its weight and thrust loads are evenly 
distributed throughout the length of the wing. If such a mounting were possible, 
it would certainly be advantageous to the conventional design. Further, assume 
that, in the proposed configuration, the cooling space between the burner and 
wing skin and the rate of flow of air in this space are such that the skin temperature 
of both wings are equal to the aerodynamic heating value for the Mach number 
being considered. This will avoid consideration of the reduced strength and 
increased boundary layer growth due to increased skin temperatui’e, since both 
wings will react the same. The wing selected for comparison was the following; 
Wing section 3% biconvex with eliptical nose 

Taper ratio 388 

Aspect ratio 3. 08 

Angle of sweep 0° 

No digedral, twist, camber, or incidence 
Values of Cjj and CX were obtained from wind tunnel tests of a wing and body 
model. Assume the proposed configuration has the following, to more closely 
match the wing used for comparison: 

Taper ratio 1 . 0 

Aspect ratio 3.08 

Angle of sweep 0° 

No digedral, twist, camber or incidence 
Wing chord 3.89 feet 



Wing span 



12. 0 feet 
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Wing Drag 

The Ct required for level flight at RI 3. 0 at 50, 000 feet is . 1305, giving 
i-i o 

CX = 6. 78° and = = 0114 for the proposed design. The closest values that 
could be found for the biconvex wing were z . 141, CX z 3. 31°, and Cq = . 0233. 
The scale factor of the model was about 20 and the effect of the change in Reynold's 
number wasn't considered. In addition, the effects of interference, helical flow 
about the configuration, and compression and expansion were not considered for 
the proposed design. The effect of considering a finite wing in the one case and 
an infinite wing in the other case is minimized since a C. was selected from 
infinite wing considerations and theCX and Cq merely transcribed for this value. 
Conservatively, then, the drag of the conventional wing is approximately twice 
as great as the drag of tlie proposed design. 

Wing Weight 

If it is assumed that the spanwise lift distribution and weight per spanwise 
foot of the two wings are approximately the same, the wings can each be broken 
at the same spanwise station and the couple loads considered. Taking the couple 
load to be 190, 000 inch pounds, a Q can be found; 190, 000 = M r dQC 
where d is the effective couple arm and C is the structural chord. The d of 
the conventional wing was taken as 1. 0 inch and the C, 50% of the local chord 
(25% to 75% MAC), to be 18. 0 inches. This gave a Q of 10550 pounds per inch. 

For the proposed design, d was taken as 5. 18 inches and C, 73% of the local 
chord (27% to 100% MAC for I'easons previously discussed), to give a Q of 1080 
pounds per inch. 

Table 4 of ref. 10 gives the required cross-section area for the upper wing 
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' surface as 6, 6296 square inches for the conventional wing. If, for the proposed 
design, it is assumed that the shock generator and burner wall carry no load 
and are equal in weight to the skin and stringers, this will give a load carrying 
cross section of 3. 33 square inches for a wing of the same weight per Lineal 
foot. This gives a Q of about 3, 000. Therefore, the proposed design of engine 
and wing the same weight as the conventional wing alone has about 2. 78 times 
the maneuver load capability of the conventional wing. 

No mention has been made of the weight of the engine and nacelle attached 
to the conventional wing, the wave drag and skin friction drag of the nacelle, 
and the interference of the nacelle and wing. Thus it can be seen that the pro- 
posed configuration is lighter by the weight of the engine and nacelle, is over 
2. 5 times as strong, provides no interference, and has much less drag (approxi- 
mately one-half the wing drag with no drag corresponding to the nacelle drag 
of the conventional configuration). The reduction in wing weight and interference 
provides an additional reduction in drag, which taken witii the previous drag 
reduction would greatly increase the range for a given amount of fuel, reduce 
the power required, and still maintain the 2. 5 strength ratio over the conventional 
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CONCLUSIONS 

From the foregoing it can be concluded that: 

1. The proposed configuration is theoretically feasible. 

2. A unique diffuser configuration can be obtained for any desired cruise 
speed. 

3. This diffuser configuration offers high total pressure recovery over a 
wide range of supersonic flight speeds, and superior design point operation. 

4. A constant diffuser throat area and Mach number can be maintained 
during operation over a wide range of flight Mach numbers. 

5. A constant burner inlet Mach number can be maintained during operation 
over a wide range of flight Mach numbers. 

6o The constant throat area and Mach number and burner inlet Mach 
number permit an optimum performance subsonic diffuser. 

7. The variable exit nozzle proposed offers improved thrust over a wide 
range of flight speeds and altitudes. 

8. The two dimensional flow considered permits engine and intake geometry 
to be varied by simple mechanical means. 

9. The integration of a ramjet power plant with a ducted lifting surface 
offers the following significant performance gains in the supersonic flight regime; 

(a) Greatly reduced drag 

(b) Stronger construction 

(c) Less weight 

10. The ramjet engine that has been proposed provides sufficient thrust to 
propel the vehicle considered for operation over a wide range of supersonic Mach 
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numbers and altitudes. 

11. The integrated system has practically constant orientation to the 
relative wind; so the variable geometry of the diffuser can be made solely a 
function of flight Mach number. 

12. These variable diffuser surfaces ai'e in part balanced by air loads; 
so that the structural requirements are not great, and the power required to 
position the surfaces is reasonable. 

13. Aeroelastic effects are greatly reduced to provide design simplification. 

14. There are a minimum number of moving parts, and none move rapidly 
enough to build up dynamic loads. 

15. These moving parts have counterparts on normal aircraft, and can be 
controlled b y single common sensing devices. 

16. It has been shown theoretically that the configuration is highly ma- 
neuverable over a wide range of flight speeds, this from the high thrust to drag 
ratio as well as the increased structural strength and lift capability. 

17. Due to the greatly reduced weight and drag, and thus thrust require- 
ments; a smaller vehicle can be built that is capable of the same range and 
payload. 

18. The design can be broken down into three areas; the supersonic 
diffuser, the engine, and the wing. Those three areas are no more difficult, 

and in most cases simpler, to design and manufacture than their normal counter- 
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Fig. 4. 

Schematic -Two Shock System 
Shock Generator and Compressor Lip Position 
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Fig. 7. 

Wing Cross Section Rear Spar Area Detail 
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Wing Cross Section at a Front Spar Post Detail 
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THESE TWO AREAS ARE MADE EQUAL TO THE FRONT 
SPAR POST FRONTAL AREA PER LINEAL FOOT OF WING 
TO COMPENSATE FOR THE POST OBSTRUCTION TO FLOW 
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Fig. 11 
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CALCULATION OF PERFORMANCE OF A TWO SHOCK SYSTEM, Mj^p 3. 0 

APPENDIX I 



The performance of a vehicle using the proposed system was calculated 
by finding the lift coefficient required for level flight of the vehicle with an 
assumed wing loading of 200 Ib./ft. ; the drag coefficient of the wing when pro- 
ducing the required Cj^; the Cj-, of the rest of the vehicle related to the wing 
chord per spanwise foot of wing; the inlet momentum term per foot of wing, 
w Vo/g, related to the wing and expressed in dimensionless coefficient form: 



w V, 



IL 



2 

M C 
o 



r- cr / ^ 

0 & / 

and the exhaust momentum term per foot of wing: 



Cf 



Tg ^ 



r/2 M^ c“ 



For level flight, then: 



C -Cp zCp_ Ci-v + C_ 

^e ^0 ^ ^Body ^wing 

since doubling the throat height doubles the wing chord; Cp , C^ . and Ci^ 

^ e ^ o ^wing 

are independent of throat height for any given configuration with constant wing 
loading. Therefore, a throat height of one inch and thus the configuration 
discussed in the text was used for calculation. This made the wing chord 
equal to three feet. 

CALCULATION OF Cp 

The station numbers used as subscripts in the following calculations 
are given in the list of symbols. The following assumptions are made: 



1 - 2 . 



P / 

Subsonic compressor efficiency, Tg / Pj 
Burner efficiency 

Specific heat of combustion products 



1.30 



. 90 



.90 



Fuel air ratio: 



f, 



max 



.071 



f . 
mm 



. 010 



Heating value of the fuel 



HV 



19, 000 BTU/lb 



No friction losses in the burner 
Isentropic flow through the exhaust nozzle 
ICAO Standard Atmosphere conditions 
Adiabatic flow through the diffuser 

The following calculations are made in Table I; The total pressure to static 
pressure ratio of the free stream flow is given by the formula: 



This ratio is given on Line 1. The ratio at the entrance to the burner will 
therefore be: 



this ratio was reduced to . 959 by the off-optimum diffuser position discussed 



Line 3, with the total temperature given on Line 4 from the formula: 





Pt 

in the text. 2 /Pq is given on Line 2. 

T 

The total temperature ratio of the adiabatic flow, T^2/Tq is given on 




1-3. 



where was taken from ICAO Tables. The increase in total temperature 



through the combustor is the product of the combustor efficiency and the theoretical 

temperature rise. This theoretical rise is a function of the inlet total temperature 

and the fuel-air ratio. In the cases underlined in Lines 5 and 6, this fuel-air ratio 

(maximum) was more than enough to accelerate the outlet air to Mg 1.0; so for 

these cases f was reduced to a value that just gave Mg 1. 0 (at the exit of the 

burner). This total temperature rise is given on Line 5 with the burner outlet 

total temperature (Line 4 plus Line 5) given on Line 6. 

For the addition of heat in a constant area duct in the absence of friction 

it is possible to obtain a relation between the total temperature ratio and the 

Mach number from a ^function or Rayleigh Line Chart. For a constant burner 

inlet Mach number the 0 function will be constant and was found to be . 0367. 

(fXMs) can then be found from the relation; 

0 (Mg) = (/) (M 2 ) ”^T 3/T-p ° ^ Line 7 

^(Mg) can then be used to re-enter the chart to find the burner exit Mach number. 

T 

The values found for the exit are given on Line 9. ^3/Trp is given on Line 7, 

^2 

and '^(Mg) values are given on Line 8. 

With an assumed height of the diffuser throat of one inch, the height of the 
burner for a constant M 2 is constant and given by the A/A* ratio for ^^of 1. 4. 

The burner height was found to be 2. 965 inches. Writing the continuity equation 
at the exit nozzle throat and the diffuser gives the ratio; ^ , 

An. m *'t 1 ‘xl ‘■'“I ,*ZlL l2(r-l) 



Mdt 


Pt| 


r Tx) 1 

'NT 


1/2 


r J 

2 


2 


Mnt 


M NT 


1 rpj 

L ^'dt - 






■"2 

“d 



Since both throats are choked, this simplifies to; 
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The first ratio was found to be . 96; and the second ratio was assumed to be l/. 90. 
found from the formula: 

^T2/P2 ^'^2/P3 

” 1 + vm| 1+ 

For constant Mg, the first ratio is a constant and was found to be . 9735. The 
second ratio was found from a chart® for a of 1. 3. These values are shown 
on Line 10. Therefore: 



P 




. 9735/Line 10 



This is sho\vn on Line 11 If one lineal foot of wing is considered: 

l/2 

12 r 1. 158 X Line 11 X (Line 7 ) ' 

p 

These values are shown on Line 12. ^3 /Pq was obtained by dividing Line 2 

by Line 11. These values are shown on Line 13. This value was then used to 

enter Thnist Coefficient vs. Area Ratio Tables to get Cj- and € for optimum 

expansion. The results are given on Lines 14 and 15. The area of the exit 

was then calculated using: 

A_ = A_C or Line 12 X Line 15 
^ 1 

The results are given on Line 16. In cases where the optimum exit area was 

limited by the thickness of the wing (those values starred), it was assumed that 

f was reduced to reduce the Mach number at the exit of the burner and thus 

reduce A to the optimum value of c for Cr . These reduced values are 
T ^max 

shown on Line 12.5. In order to simplify calculations it was assumed tliat A , 







r ' 





.f: 




* J 






i 



i 
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the correction factor for the divergence angle of the exhaust nozzle was unity; 
and that even though f was reduced greatly, the total pressure and total temperature 
loss remained the same. It was felt that this was very conservative. In the 
range 3, 0 and greater this f value required for optimum € got below the 
assumed minimum fuel flow. Values for f of . 010 were therefore recalculated 
in the short table after Line 16, and are shown in the table starting with Line 17. 
The values of Cj found for this condition are shown in the table. It should be 
noted that they range from 99.6% to 97. 7% of the previous optimum C^. 

The exhaust thrust divided by the static pressure was found from; 

TT C' P 

■^e f 12 A't , Tq Line 12 X Line 13 X Line 14 

— = 12 —^ »■— 12 

(2 

The results are given on Line 17. This was converted to on Line 18. 



'Fq was computed from the equation: 



Cp _ vV Vq 
o 



g 



/ 



r 



w was found by finding a constant: 

w-Vel 






y+1 



A - 



^2 St^ ■ 

where; 0. = (T^ 

to 'A; 



M, 



(1 + 4 -^ «2 > ^**”'^* 



2/T„) L. ) 



6 



‘2 - > ‘Palt/Pg > 

Values of (5j.^ are given on Line 19, values of are given on Line 20; and 



values of w are given on Line 21. C„ will be the same at all altitudes; so only 

^o 

the values of V /g for sea level are given on Line 22. C„ values are given 
^ ^o 

on Line 23. The engine coefficient of thrust is therefore: 



C 



F = 




- C 



F 



o 



or Line 18 - Line 23. 




» 


















M 

oJ . ’ 



1 




t 



1^0 
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TMs is shown on Line 24 and in Fig. 10, 

CALCULATION OF C^ AND Cj^ 

The Cj^ required for level flight was found by the formula; 

Ct = ^ where VV/S = 200 lb. /ft. ^ 

^ ^/2 P ^ m 2 

o o 

These values are shown on Line 25. 

Showing the thrust to be well over the drag for one altitude demonstrates 
the ability of the vehicle to accelerate to cruise speed. Using a throat height 
of two inches makes all the cross-section dimensions twice as big as those 
discussed in the text. This makes the actual chord: 

lo625 ft. (generator) + 4.375 ft. (aft wing) = 6 ft. 
generator chord = .271 total chord 
aft wing chord = .729 total chord 

The height of the throat is limited by the length of generator that can be 
controlled and held in place; but 1. 625 ft. seemed very reasonable. The Cj^ and 
Cj) of the shock generator were computed from values of 0^ given on Fig. 9, 
using Linearized Thin Airfoil Theory, The results are given on Fig. 9 and in Lines 
26 through 30. The body shape discussed in the text with of 1. 38 ft. was 

assumed. The Cj^ values related to the wing chord are given on Line 31, \vith 
the total Cjj given on Line 32, and in Fig. 10. Lines 26 through 30 assumed that 
the vehicle could be rolled 180° (turned upside-down) for cruise in order to keep 
the angle of attack nearly constant at zero degrees. The short table following 
Line 32 was used to investigate the effect of rolling the vehicle at 3. 0 
and greater speeds. It can be seen in botli cases that Cp is well above C^. 
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In order to simplify calculations, the lift of the body and tail surfaces, and 
the relatively large portion of the wing submerged in the fuselage were ignored. 
This would be offset somewhat by the non-linear effects of interference of wing 
and wing, and wing and body that were also ignored, but the simplification will 
still be conservative. The calculations were further simplified by making a 
pessimistic estimate of the lumped drag effects of body, tail, and submerged 
wing. 
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TAriLil, I~1 f 

CALCULATION OF OURFORL^ANCU OF A TUG SHOCK SYSTEIvi, M^p 3.0 
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2.6 


3.0 
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Alt. 
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1 D • ^ (5 


10.7 


19.95 


36.7 
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1 


^T2/Eo 


3.2 


I 4.65 


9.3 


17.01 


29.9 
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1 78.4 


— 


2 
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2.80 


3.45 
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.1092 


1 .0931 


.0814 |0 K 




I 

' o 


.2172 1 


.2172 


.1637 


^ i 

.141 1 


1 

.1225 


.1040 


.090 1 

1 


I 20 K 


8 


\ 

1? 

1 


. 2 I 72 I 


.2172 


.2172 


.1589! 


.1369 


1 .1156 


< 
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SHORT TABLE 1-2 



0 


K 
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3.0 


3.5 


3.0 


3.5 


3.0 


3.5 


He 
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■^3 
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1 1.43 


1.35 


1.49 
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.2385 


.251 
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“3 
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. 954 


.966 


.962 
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1 

Prj / 

f^/ • +^^2 


1.007 


1.006 


1.01 


1.005 
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1.009 




1.34 


1.328 


1.4 


1.354 


1.43 
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-^NTl 














■P 


29.7 


50.9 


29.6 


50 .6 


29.55 
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2.985 


3.01 


2.86 


2.95 


2.8 


2.S9 


e 
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1.463 


1.441 


• 

CD 


1.439 


1.476 




4.79 


8.35 

/ 


4.93 


&.46 j 


5.07 


e .62 


5’«/Po j 


. 2535 
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.335 





Mo 
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3.0 
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• 
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H 
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1.05 


J^.195 


3.375 


30e 


above 


0 K 


1 

j 


F /P 

9 ' 0 


.567 i 
1 


i ,950 


2.190 


3.34 


see 


above 


20K 


! 17 
\ 




1 

.567 


\ . 960 
) 


2.045 


3 .32 


see 


above 


36K 


1 

\ 




.120 


.1635 


.216 


.237 


see 


above 


0 K 






.120 


.1495 


.2155 


.235 


see 


above 


20K 


18 




.120 


.1495 


.211 


. 2335 


see 


above 


36K 






3.2 


4.65 


9.3 


17.01 


29.9 


51.2 


0 K 






1 .47 


2 . 135 


4 . 27 


7 , 81 


13.62 


OS - 
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^2. 
















19 




.719 


1.042 1 


I 2 .O 65 


3.82 


6.72 


11.47 


36K 






.368 


.535 


1.068 


1.955 


3.438 


5.88 


5 OK 


1 
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I-ll. 



Mo 1 


' 1.5 ^ 


1.75 


2.2 ' 

1 


2. '6 


3.0 ! 

1 


3.5 

1 


Ait . 


Line 




1 . 207 


1.27 


1.402 1 


1.5o3 


1.672 


1.858 


OK 




i.il9 


1.17 9 


1.304 


1 .425 
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20K 
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1 
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1 0 K 
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0 
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36.3 


56,4 


20K 


21 


2.64 


•-7 0 r? 
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29.35 


36K 
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3.62 


6.14 
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15.07 

1 


50K 
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51.9 


60,6 


76.1 

^ 1 


90 i 103.9 ! 

L 


' 121.7 


0 K 


22 




.0571 


.0676 
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.0972 I 


[ i 

.1373 j 


1 ! 

j .1925 j 


.2560 


- 


23 




.0629 


• 095 5 
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^ 1 

|.0997 i 


1 1 
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.0685 


0 K 
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CALCULATION OF PERFORMANCE OF A THREE SHOCK SYSTEM 4.0 

APPENDIX n 

For this system the throat height was assumed to be one inch, giving the 

schematic shown in Fig. 11 with values as shown in Fig. 12. This gave a 

diffuser chord of 27 inches, an aft wing chord of 27 inches, and a total chord of 

2 

54 inches or 4.5 feet. The wing loading was again assumed to be 200 lb. /ft. . 

The engine span was taken as 12 ft. with the same size body, making the actual 
wing span 14. 76 ft. and the aspect ratio 3. 28. It was again assumed that the 
missile could be rolled at cruise speed. The effect of restricting this roll was 
not considered since it was shown to be small in Appendix I. 

The calculations were made in the same manner as Appendix I, using the 
same assumptions. For speeds of M^ 3. 5 and below. Lines 1 and 3 through 12 

p 

of Appendix I were used. Line 2, the T£/p for this diffuser as shown in 
Fig. 13, with the ^ as given in that figure, is given in Table U-1, As before, 
an altitude of 20, 000 ft. was used for computations below cruise speed. These 
computations were made in Table U-1, with line numbers corresponding to the 
line numbers of Appendix I. 

Mq 4. 0 values were computed using minimum fuel flow. These values 

are given in Short Table U-2, which was extended in order to show the cruise 

condition at various altitudes. The Cj^ and Cj^ of the generator (Lines 32 and 

33) are shown in Fig. 13, and the total Cj^ and Cp (Lines 39 and 24) are shown 

in Fig. 14. It should be noted that due to the relatively thicker wing, when compared 

to the burner height; a larger C was available which gave a Cf at M^ 4. 0 of 

^ax 



n-2. 



1. 596. This was limited very slightly by wing size, and was found to be 
99. 7% of Cf optimum. As in Appendix I it can be seen that C is well above 

I X 



s ^ ^=1 
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TABLE II-l 



CALCULATIOl'- OF PE-^FORuIANCt OF a JdOCK SYSTEM 11 4.0 
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LOSSES CONSIDERED IN THE CALCULATION OF THE 

TOTAL PRESSURE RECOVERY 
APPENDIX m 

The losses considered in the computation of the total pressure recovery 

were; oblique shock losses in the convergent section, normal shock at the throat, 

boundary layer build-up in the convergent section, skin friction in the throat, 

and subsonic diffuser losses. The oblique shock losses were minimized by the 

P^ 

method of Oswatitsch, and taken into account as ^i/p hi Figs, 3 and 12. 

”^0 

The normal shock losses at the throat were minimized by the isentropic section 
of the supersonic diffuser, since this section reduces the flow Mach number so 
that it approaches Mach 1. 0 just prior to the normal shock. 

The boundary layer build-up in the convergent section is controlled in two 
ways. First, the boundary layer bleed in Fig. 6 of the original proposal will 
reduce losses on the upper surface of the convergent section. Second, wind 
tunnel experimentation has shown that the boundary layer on the lower surface 
can be minimized by holes drilled at an angle in the lower surface. These 
holes are positioned so that the number exposed to the flow will be increased 
in proportion to the tendency for the boundary layer to grow as the lower lip is 
extended. It has been demonstrated that this procedure will reduce these losses 
to negligible values. It should also be noted that the boundary layer bled in 
Fig. G is used for burner area cooling. 

Skin friction losses in the throat are incurred throughout the length of 
the constant area throat required to reduce the flow to subsonic speed. This 
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m-2. 



loss can be minimized by minimizing the required throat length. This length 
is made up of two values; a "lost length" of approximately one hydraulic diameter 
in front of the beginning of the normal shock and a length equal to the "length" of 
the normal shock., The latter length approaches zero as the Mach number before 
the normal shock approaches 1. 0. This theory has been partially confirmed 
in a wind tunnel and an average value of skin friction has been found to be . 0010 to 
. 0015. Thus the constant area throat required can be reduced to one hydraulic 
diameter and the friction loss reduced to a negligible value. 

The subsonic diffuser losses can be minimized by the use of the optimum 
expansion angle in the diffuser. This was done and a pessimistic efficiency of 
. 90 was assigned. 

This tiieory of a required constant area throat length that is related to the 

Mach number before the normal shock was then applied to the oblique plus 

normal shock diffusers. The Mach number before the normal shock of the one, 

two, and three oblique shock systems was taken from ref. 6. If it is assumed 

that the throat length is related to this Mach number, M^, as shown in Fig. III-l 

and that the throat height is one inch, with an assumed coefficient of skin 

friction, Cgp of .0015; then a total pressure loss due to skin friction in the 

throat can be computed as shown on Fig. ni-l. This loss, , was 

^o 

applied to the total pressure recovery curve in ref. 6, and the result shown in 
Fig. 3o This loss was found to vary from .23% for the proposed configuration 
to almost 2. 85% for a one oblique plus normal shock diffuser operating at Mach 
4.0. 



The above loss would seem to be negligible; but from the following, it can 
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be seen that this is not the case: 

Consider the loss as a perturbation on the thrust equation for an ideal 
expansion: 



F = 



P = P 
^e o 



where: Fq is the inlet momentum, assumed constant (the free stream conditions 

remain constant) 

Fg = Cf P^ is the exit momentum 
Cf is the thrust coefficient of exhaust expansion 
is the area of the exhaust nozzle throat 



PTc - (FTi/p^^ ) (Ft 2/ P^^ ) ( 



P Pt 

T 2 /p^ was assumed to be . 90, and ■‘■S/P.p^ a constant 

Thus it can be seen that F is proportional to . 90 A 



Assume a reasonable relationship between F^ and F^ for the cruise condition, 



say: 

therefore: 



Fo = 5/6 Fe 



F z 1/6 Fg 
■ ^p- = (6) (• 90) A^Ti/p 

Ti/p^^) Tg 

Thus the loss in net thrust for the above range of A Tl/P„ is actually 1. 24% 



for the proposed configuration to 15. 4% for the oblique plus normal shock diffuser. 
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